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Abstract 

We develop a systematic coarse-graining procedure wtiicli establislies tlie connection between 
models of mixtures of immiscible fluids at different length and time scales. We start from the 
Cahn-Hilliard model of spinodal decomposition in a binary fluid mixture under flow from which we 
derive the coarse-grained description. The crucial step in this procedure is to identify the relevant 
coarse-grained variables and find the appropriate mapping which expresses them in terms of the 
more microscopic variables. In order to capture the physics of the Doi-Ohta level, we introduce the 
interfacial width as an additional variable at that level. In this way, we account for the stretching 
of the interface under flow and derive analytically the convective behavior of the relevant coarse- 
grained variables, which in the long wavelength limit recovers the familiar phenomenological Doi- 
Ohta model. In addition, we obtain the expression for the interfacial tension in terms of the 
Cahn-Hilliard parameters as a direct result of the developed coarse-graining procedure. Finally, 
by analyzing the numerical results obtained from the simulations on the Cahn-Hilliard level, we 
discuss that dissipative processes at the Doi-Ohta level are of the same origin as in the Cahn- 
Hilliard model. The way to estimate the interface relaxation times of the Doi-Ohta model from 
the underlying morphology dynamics simulated at the Cahn-Hilliard level is established. 
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I. INTRODUCTION 



The problem of coarse graining in terms of bridging time and length scales between 
microscopic and macroscopic levels of description is a crucial issue in the physics of com- 
plex fluids, hke polymer melts, colloids, liquid crystals and emulsions, out of equilibrium. 
The wide span of length and time scales is particularly evident in these systems, due to 
their internal structure leading to additional mesoscopic levels of description, intermediate 
between microscopic and macroscopic [1]. Many of the practical applications of these flu- 
ids crucially depend on the evolution of their complex, multiphase morphologies developed 
through equilibrium self-assembling or during non-equilibrium processing of these systems. 
Some examples of such applications are food processing, membrane technology, encapsula- 
tion systems, drug delivery, coating, production of paint and cosmetics etc jl-4|. Therefore, 
the goal of connecting different levels of description of complex fluids is of great importance. 
In the present work, this problem is approached by considering the phase separation of 
binary fluid mixtures subjected to a shear flow. 

When a binary (AB) fluid mixture is quenched from a high temperature homogeneous 
phase to an unstable region below the coexistence curve, it becomes unstable with respect 
to long- wavelength fluctuations in the composition and starts to phase separate jsj. The 
interface between the A- and B-rich domains is initially very diffuse, but sharpens with time. 
In the later stages of phase separation, local equilibrium is achieved within each domain and 
the effects of interfacial energy become important. The domain pattern further coarsens in 
time, while the dynamics is governed by the minimization of the excess interfacial energy of 
;he system. The kinetics of this nonequilibrium phenomena is an area of extensive research 
5| . Especially the effects of a shear flow on the deformation and kinetics of domain growth, 
and the corresponding rheological behavior, are challenging and technologically important 
topics jo]. The complexity of the problem lies in the presence of a complex interface and 
its motions due to coagulation, rupture and deformation of domains, which significantly 
influences the macroscopic properties of a mixture. Although numerous experimental 
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numerical 12N24| and analytical 
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17 



24 



28| studies have been done in order to 



understand the flow effects on the morphology and the dynamics in this system, many 
questions still remain open, like the possibility of achieving a nonequilibrium steady state 
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-|24|. 
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Due to the complex morphology of this system, even in the case of two simple Newto- 
nian fluids, one can identify several length and time scales requiring different theoretical 
approaches and simulation methods. The standard mesoscopic model for understanding the 



dynamics of spinodal decomposition is based on the formulation of Cahn and Hilliard |29 |. 
It describes the kinetics of the system in terms of the convective-diffusion equation for the 
order parameter - composition c, and the Navier-Stokes equation for the fluid velocity v 



(so-called 'model H', see e.g. 



3Q|): 



dc 

+ v(Vc) = MVVc, (1.1a) 

ot 

P + (v ■ V) = r/VV - Vp + /icVc, (1.1b) 

where M is the mobility coefficient, p the pressure, p the mass density, and rj the viscosity. 
The chemical potential pc is given as pc = SF [c] / 6c, with the free energy functional of the 
Ginzburg-Landau form F [c] = J d^r {-{a/2y + (6/4)c^ + (k/2) \Vcf), (a, 6 > 0), where 
the square brackets emphasize the occurence of functional integrations. During the late 
stages of phase separation, the time evolution can be described by the change of the size 
and shapes of the domains. Due to the dynamical scaling hypothesis, the system is described 
by a characteristic length scale L{t), which is related to the average domain size. By using 
dimensional arguments, one can estimate the size of the contributions from the specific 
ierms in the equations (11. ip . which leads to the distinction of the three growth regimes 



E 



32|: diffusive L(t) ~ (Mat)^^^, viscous hydrodynamics L{t) ~ o-t/rj, and inertial 

1 /3 

hydrodynamics L{t) ~ [ot^ / p) , where cr is the interfacial tension. Under shear flow 
the situation becomes more complex since the domains elongate along the flow direction 
and the system becomes anisotropic and can not be described by a single length scale (see 
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20| and references therein). 
On the other hand, for understanding the rheological properties of multiphase systems, 
like polymer blends, under shear flow, Doi and Ohta suggested that the average information 
about the interface between the two phases is sufficient. Hence, their phenomeno logical 
model 28| focus on the interface, which is considered as a zero width surface embedded in 
the fluid. The presence of the interface is then represented through two state variables that 
describe the interfacial area per unit volume, Q, and its anisotropy, q, in a given flow fleld. 
They are deflned in terms of the interface orientational distribution function / (n), where n 
is a unit vector normal to the interface. This is a probability density of flnding the amount 



of interface with the orientation n, so that the integral over all orientations gives the total 
amount of interface per unit volume. Therefore, for the state variables Q and q one can 
write the following definitions 



Q 
q 



nn - ^1 j / (n) (fn = nn - ^Ql, 



(1.2a) 
;i.2b) 



where 1 is the unit tensor, and j (Pn denotes an integration over the unit sphere. Introduc- 
tion of the distribution function and further averaging over all possible orientations imply 
the loss of information of the detailed morphology of the system. This coarse-grained model 
keeps only the information about the average amount of interfacial area per unit volume and 
its orientation. However, the detailed information about the morphology, i.e., the explicit 
interfacial position and orientation, is lost. 

The time evolution of the new configurational variables Q and q is determined by two 
factors: the external flow field which orients and enlarges the interface, and the interfacial 
tension which has the opposite effect and governs the relaxation of the interface. Therefore, 
the time evolution can be separated into two parts: 

dQ_ 
dt 
9q 

convection i 

The convective part of the time evolution can be found by considering the convection be- 
havior of the unit vector n normal to the interface under affine deformations and volume 
preservation assumptions. Then, the convection of the state variables is derived using the 
above definitions fll.2D 



dQ 


dQ 


dt 


dt 


d(\ 


9q 


'dt ^ 


'dt 



+ 



convection 



relaxation 



relaxation 



(1.3a) 
(1.3b) 



dQ_ 
dt 
9q 
'dt 



convection 



convection 



-V-(Qv)-(Vv)^:q+^Ql:(Vv)^, 

- (Vv) ■ q - q ■ (Vv)^ _ ig^ _ V ■ Vq + ^1 (Vv) : q 



;i.4a) 



+ 



Tr (Vv) + nnnn : (Vv)""" . 



;i.4b) 



where 7 = Vv+ (Vv)^ is the symmetrized velocity gradient tensor [33|]. In the equation for 
the interfacial shape ( ]1.4b|) . the fourth moment of the interface normal n appears, which must 
be expressed in terms of the state variables through an appropriate closure approximation 



in order to obtain a self-contained set of time evolution equations. Doi and Ohta postulated 
the following closure approximation 



nnnn : (Vv)^ = — nn nn : (Vv 



:i.5) 



and verified its accuracy. Using definition (ll.2b|) . the convective time evolution of the 
anisotropy q takes the form 



'dt 



convection 



(Vv) ■ q - q ■ (Vv)^ _ ig^ _ V ■ Vq + ^1 (Vv) : q 



+ 



q 



-Qi 

9^ 



Tr (Vv) + ^ (q + Iq^J (q + ^Ql ) : (Vv) . (1.6) 



As for the relaxation in the Doi-Ohta model, the relaxation of the state variables can be 
written in a generalized form 

= - -^Q, (1.7a) 

relaxation ^do,l 

^ -q. (1.7b) 



dQ 
dt 
dq 
'dt 



relaxation ^do,2 

In order to specify the empirical interface relaxation times rdo,i and rdo,2, assumptions on 
the appropriate interface relaxation mechanisms must be made. Finally, the full set of time 
evolution equations for the configurational variables Q and q, together with the equations for 
;he hydrodynamic variables p, g and e has been formulated within the GENERIC formalism 



3^ . which for the stress tensor gives 



cr{ci--Ql] -pi 



where r is the viscous stress tensor. 

Although very crude, the ability of the Doi-Ohta model to connect the morphology of the 
system to its macroscopic behavior, through the above given set of rheological constitutive 
equations, is very appealing for different theoretical and technological purposes. The model 
and its variants have been used to analyze multiphase flow of polymers, as well as simple 
fluids undergoing spinodal decomposition 0, [s^^O]. However, the need for explicit assump- 
tions of the interfacial relaxation mechanisms and the appropriate relaxation times, which 
depend in an unknown manner on the underlying morphology, makes it difficult for wide 
and easy use. Therefore, in this paper we concentrate on relating the macroscopic Doi-Ohta 
model to the underlying morphology, described by the more detailed Cahn-Hilliard model. 



The goal is to establish a systematic, thermodynamically guided method which determines 
the coarse-grained Doi-Ohta model of a phase separating binary fluid from the more detailed 
Cahn-Hilliard model, in the region where the considered system can be described with both 
approaches. To connect these two different levels of description, the General Equation for the 
Non-Equilibrium Reversible-Irreversible Coupling (GENERIC) framework is used 4lN43|. 
This approach allows to specify the mathematical structure which should be left invariant 
in the coarse-graining procedure and hence assures thermodynamically consistent results. 

In the following, we flrst describe the GENERIC structure and formulate the Cahn- 
Hilliard model within this framework. The coarse-graining procedure which derives the more 
macroscopic Doi-Ohta model is then presented. In particular, we present some mathematical 
challenges which occur within the coarse-graining procedure, as well as the efficient way to 
extract the interface relaxation times of the Doi-Ohta level from the underlying morphology. 



II. GENERIC FORMALISM 



The main points of the GENERIC framework of nonequilibrium thermodynamics |4l|-|43| 
can be summarized in the following way. In analogy to equilibrium thermodynamics, one 
needs to choose a complete set of variables x, which describes the situation of interest to the 
desired detail. For a system out of equilibrium the time evolution of all the relevant state 
variables can be divided into the reversible and the irreversible contributions 

X x|i-Qv -|- xjiirev • 

These two contributions are obtained by considering separately the energy E and the entropy 
S — the two generators of reversible and irreversible dynamics, respectively. This can be 
explained by using the analogy with the classical Hamiltonian mechanics as following. Using 
the classical Hamiltonian mechanics, the reversible contribution of the time evolution, x|rev, 
is related to the energy gradient by way of a Poisson operator L, so that it is xlrev = 
L ■ {6E/6x.). The Poisson bracket associated to L, given by {A,B} = (5A/5x, L ■ SB/S:k.) 
with appropriate scalar product (.,.), must be antisymmetric and satisfy Lebniz' rule and 
the Jacobi identity, which are both features that capture the nature of reversibility. The 
motivation for formulating the irreversible part of the dynamics comes from the reversible 
part. For the reversible dynamics, the energy plays a distinct role — it is a conserved quantity 
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for closed systems and it drives the reversible dynamics. Therefore, for the irreversible 
dynamics the entropy presents an important quantity — it does not decrease for closed 
systems. For the GENERIC formulation, it is then assumed that the irreversible contribution 
to the time evolution, xlirrcvj is driven by the entropy gradient, i.e., that it is of the form 
x|irrev = M ■ (5S'/5x), with M a generalized friction matrix. This friction matrix contains 
transport coefficients and relaxation times associated to the corresponding dissipative effects. 
It is required for M to be (Onsager-Casimir) symmetric and the condition that it is positive 
semi-definite ensures that > is fulfilled. 

Prom the above illustration, one concludes that the time evolution of x can be expressed 
in terms of four "building blocks" E, S, L and M as 

X = L(x) • + M(x) • . (2.2) 

The two different contributions to the time evolution, reversible and irreversible, are not 
independent. They are related by the two complementary degeneracy requirements 

L(x).^f^ = 0. M(x).^f^=0. (2.3) 

OX ox 

The first requirement expresses the reversible nature of the L contribution to the dynamics, 
demonstrating the fact that the reversible dynamics captured in L does not affect the entropy 
functional. The second one expresses the conservation of the total energy of an isolated 
system by the irreversible contribution to the system dynamics captured in M. 

The presented GENERIC form of the time evolution equations for the state variables x 
represents a mathematical structure which guarantees that the chosen model is thermody- 
namically consistent. Moreover, by concentrating on the thermodynamics building blocks E, 
S, L and M rather than on the time evolution equations in systematic coarse-graining pro- 
cedures, the obtained coarse-grained equations are guaranteed to be also thermodynamically 
admissible. 



A. GENERIC formulation of the Cahn-Hilliard model 

In order to develop a systematic coarse-graining procedure which relates the more detailed 
Cahn-HiUiard level to the more macroscopic Doi-Ohta level, we must formulate the first one 
in the GENERIC structure. As a first step, we need to identify the list of independent 



7 



variables x which fully describe the considered thermodynamic system. In addition to the 
hydrodynamic fields, the mass density p, the momentum density g, and the internal energy 
density e, we use the composition c - the mass fraction of one component - in order to 
describe the configuration of the system. 

Taking into account the Cahn-Hilliard free energy of a phase separating binary system, we 
assume the total energy and entropy of the system, in terms of the variables x = (p, g, e, c), 
to be 



E«(x) 
5«(x) 



V2 p 

d^r s (p, e, c) 



-ke I Vc| 



(2.4a) 
(2.4b) 



where V is the total volume of the system. For simplicity, the gradient-squared term, which 
describes the contribution due to the presence of the interface, is included only in the en ergy 



functional fl2.4ap . although it can be separated between both, energy and entropy 4J, |45 |. 
The actual functional form of the entropy density s (p, e, c) is not needed for the derivation 
of the Cahn-Hilliard time evolution equations in the GENERIC form. However, we use 
the assumption of local equilibrium, i.e., s (p, e, c) has the same form as an equilibrium 
system at the corresponding state point. In the numerical simulations reported in this 
paper, we have used the popular c^-structure for the uniform part of the entropy density. 



i.e., s(p,e,c) = Sq (p^ 



free energy density 
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ac^ 1 2 -|- hc^ lA, which presents the entropic contribution to the bulk 
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19 



20 



22 



24j |. Other possible functional forms include the one 



for the van der Waals systems, for example Eq. 4 in [4j|. The functional derivatives of the 
above energy and entropy functionals with respect to the independent variables x are given 
with 

T 



5x 



— v^v, 1 

2 ' ' 

ds (p, e, c) 



dp 



,0, 



keV^c 

1_ ds (p, e, c) 
T' 



dc 



(2.5a) 
(2.5b) 



where we omitted the position dependence for simplicity. 

The Poisson operator L*^^) determines the reversible contributions to the full evolution 
equations of the variables x, and hence models the convective behavior. For our particular 
choice of variables and in view of the functional derivatives (12. 5 p one arrives at the following 
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expression for the Poisson operator 

/ 



L(i)(r,r') 



P (r') ^ 








L 



(1) 

99 

-(1) 
"^9 



L 



(1) 

9^ 





\ 

(Vc(r))(5 




(Vc(r))(5 







(2.6a) 



with 



r, r 



9(5 95 
ot' ar 

OT or 

e(r) — +p(r) — , 



dr' 



(2.6b) 
(2.6c) 
(2.6d) 



where 6 = 6 {r — r') is Dirac's 5-function. In addition to the usual hydrodynamic part 
of the Poisson matrix for the variables p, g, e, element L^cg^ dictates the convection of 
the configurational variable c, which is a scalar, so the entry in the Poisson operator is 
as expected by The antisymmetry of the Poisson matrix fl2.6p and the degeneracy 

requirement are satisfied by construction, and the Jacobi identity holds. At this point we 
mention that, in general, the symbol "■" in (12. 2p implies not only summation over discrete 
indices. If field variables are involved the operators L and M are written in terms of two 
space arguments (r, r'), and an integration over r' must be performed when multiplied with 
a fu„c«o„ of r' from the right. However. i,r the case of the field equations being loeal. one 
can express L and M in terms of a single variable r only [41|, and no integration is implied 
when these operators are multiplied from the right. Such single variable notation is used 
below for the friction matrix M*^^-'. 

The irreversible effects in the Cahn-Hilliard model for a binary fluid under flow arise due 
to viscosity and heat conduction, as well as diffusion. Thus, to incorporate them through the 
friction matrix M*^^\ one can write it as the sum of two matrices M*^-*^-' = M*^-'^^'^^'^ + M*^-'^'''*^'^, 
where M*^^-'''^^'^ is the usual hydrodynamic friction matrix, and M*-^-'''^'^ contains all the 
transport coefficients related to diffusion of mass 4l|, |46|. The hydrodynamic part of the 
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friction matrix is given in the r-notation as 



/o 








o\ 

















^l).hyd 













0/ 



(2.7a) 



with 



- (Vr^TV + IV ■ riTVf - VkTV, 
V ■ T]Tj + V— tr7, 

-r]Tj ■ V - — tr7V, 

f 7 : 7 + ^ (tr7)^ - V ■ A^TV, 



(2.7b) 
(2.7c) 

(2.7d) 

(2.7e) 



with all derivative operators acting on everything to the right, i.e., also on functions multi- 
plied to the right side of the operator M, and with with the transport coefficient k, being 
a combination of the dilatational viscosity k, and the viscosity r], k = k — |?7, and A'' the 
thermal conductivity 



If M (r) is a position dependent mobility coefficient, the diffusion part of the friction 
matrix in the r-notation is 

/O \ 





M(i)'dif (r) 



mL'^''^'' MiV'"'' 
\0 M^^''"'' Mg^''''' J 



(2.8a) 



where the cc-element is describing the diffusion of c, so that 

M(l),dif ^j.) _ _y . ^j^^^ ^2.8b) 

and the other three elements are obtained from the symmetry and the degeneracy conditions 

M(i)'dif (j.) = _v . MTVke (V'c) , (2.8c) 
-ke (V^c) V ■ MTV, (2.8d) 
-ke (V^c) V ■ MTVke (V^c) . (2.8e) 



M(i)'dif (r) 



Here again are all variables functions of r, and all derivative operators V = ^ acting on 
everything to the right, except when they are in brackets and bounded to act on c. The 
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symmetry and the degeneracy requirement for the matrix M*^^ ^''^'^ are satisfied by construc- 
tion. 

Finally, by inserting the above building blocks E, S, L and M into the GENERIC 
equation (12. 2p . the full set of time evolution equations takes the form 



dg 

dt 



-V ■ (vg) - V ■ (n + r) , 



•9e ^ , . 



pV-v-r : (Vv)^- V-j 



dc 
di 



+ KE (V^c) V ■ 

-V ■ (Vc) + V 



MTV 
MTV 



ds{p,e,c) ke^2 



dc 
ds (p, e, c 



dc 



(2.9a) 
(2.9b) 

(2.9c) 
(2.9d) 



where 11 = (p — |Vc|^) 1 + ke (Vc) (Vc) is the pressure tensor, t = ri'y — k (V ■ v) 1 is 
the viscous stress tensor, and j*^ = — VT represents the conductive fiow of internal energy. 
For an isothermal, incompressible fiow, these time evolution equations take the form of the 
standard Cahn-Hilliard model. Similar set of hydrodynamic equations for phase separating 



fiuid mixtures has been derived from an underlying microscopic dynamics 



46|. 



III. COARSE GRAINING FROM CAHN-HILLIARD TO DOI-OHTA LEVEL 



A specific feature of the GENERIC formalism is that it is applicable at different levels 
of description associated with different length and time scales. Each level of description is 
described by an appropriate set of state variables and building blocks. The mathematical 
structure of GENERIC offers a possibility to perform coarse graining by focusing not on the 
time evolution equations, but on the building blocks E, S, L and M, see The transition 
between a more detailed (Cahn-Hilliard) level 1 to a less detailed (Doi-Ohta) level 2, which 
involve the derivation of the building blocks from level 1 to level 2, can then be performed 
by systematic procedures 48|| . The possibility to perform coarse-graining by focusing on the 
basic building blocks guarantees that the thermodynamic structure of the problem will be 
preserved. 

A systematic coarse-graining procedure is based on the idea of the separation of time 
scales. In this way, it is assumed that there exists a division of the "fast" and "slow" 
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degrees of freedom at the more microscopic level 1 according to their relaxation times. By 
means of the projection-operator technique, one can then eliminate these "fast" degrees 
of freedom from the time evolution equations for the "slow" variables. The latter are then 
associated to the macroscopic variables at level 2. The crucial steps in this general procedure 
are to identify the proper mapping of the variables of one level to another, n(x) : x — > y, 
which average in a non-equilibrium ensemble py (x) is the new coarse-grained variable y 



n(x) ) = / Vxpy[x]U{x) . (3.1) 



A. Mappings and ensemble 

For coarse graining from Cahn-Hilliard level, the relevant set of variables at the coarse- 
grained level y = (p, g, e, Q, q) is motivated by the phenomenological Doi-Ohta model 
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33| . We assume that the hydrodynamic fields p, g, and e are smooth on the more 
microscopic Cahn-Hilliard length scale. Then the mappings Up, Ilg, and 11^, simply pick 
out the hydrodynamic fields, while we need to determine the relationship of mappings IIq 
and IIq to the underlying configuration. We do this by following the physical meaning rather 
than the exact definitions (11. 2p of the Doi-Ohta variables. This is because, first, the interface 
orientational distribution function / (n) is not given at the Cahn-Hilliard level, and, second, 
due to the difference in modeling of the interface. While in the Doi-Ohta model the interface 
us assumed to be sharp, in the diffuse-interface theories, like Cahn-Hilliard one, the interface 
is given through continuous variations of the composition c. Therefore, for determining the 
average interfacial area per unit volume Q and its orientation q one must use the appropri- 
ate combination of the gradients of the composition Vc(r), and perform ensemble averaging 
which incorporates "smoothing" over a certain volume. We introduce a smoothing function 
x(r — r'), which averages the observable over a certain volume v{r) around position r, and 
satisfies the normalization condition jyd^r' x(r — r') = 1. To obtain the variables which 
would correspond to the Doi-Ohta averaged interfacial area and its orientation, volume v{r) 
must comprise many droplets for good statistics. That means that the length scale of the 
smoothing volume v{r) - smoothing length a - must satisfy 

e < L{t) < a < A, (3.2) 
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where ^ is a length of the size of the interfacial width, L{t) is the growing characteristic 
domain size, and A is the size of the system. The first part of the inequahty (13.21) . C, ^ L{t), 
denotes an important fact that we perform the coarse-graining procedure only during the 
late stages of phase separation, i.e., when the well defined domains are formed and the 
coarsening can be described by the growth laws given in Section [H 

Although one can express the coarse-grained variables Q and q using gradient of the com- 
position Vc(r) and average the expression as discussed above, there are still some questions 
which are not solved. As already mentioned, contrary to the phenomenological Doi-Ohta 
model, in the coarse-grained model, the interface not only has a finite width, but it does 
not behave like a zero-width mathematical surface. Rather it deforms under the applied 
flow and its stretching must be captured by the new variables, which is not the case in the 
phenomenological Doi-Ohta model. There are different ways to solve these problems. One 
way would be that instead of using the Cahn-Hilliard Poisson operator L^^\ Eq. (12.61) . one 
needs to use its modification with the appropriate constraint which would account for the 
stretching of the interface. Another possibility would be to make certain modifications to 
the original Cahn-Hilliard model in such a way to make the interfacial width, as well as in- 
terfacial tension, fixed. However, maybe the most straightforward way to solve the presented 
problems is to account for the finite interfacial width and its stretching by introducing an 
additional variable into our coarse-grained Doi-Ohta model. This additional variable would 
be the average interfacial width /. This means that for the configurational variables of the 
coarse-grained Doi-Ohta model we use new variables {P, p, /}, which are obtained as ensem- 
ble averages of the suitable mappings Hp, Hp and H;. The original Doi-Ohta model and 
the appropriate convective behavior of its variables are obtained by transformation of the 
configurational variables {P(r), p(r), /(r)} to {(^(r), q(r), /(r)} as 

g(r) = /(r)P(r), q(r) = /(r)p(r). (3.3) 

We choose the appropriate mappings for the new variables, having in mind that the 
variables Q and q, obtained from {P, p, /} as (13.31) . should represent the average amount of 
interface per unit volume and its orientation, both of dimension m~^. Then, for suitable 



13 



mappings for the new variables we propose 

np[c](r) = /" rfV |Vc(r')|^x(r-r), (3-4a) 
Jv 

np[c](r) = ^rfV (^(Vc(r'))(Vc(r'))-i|Vc(r')|'l)x(r-r'), (3.4b) 

! rfV |Vc(r')|x(r-r') 

n;[c](r) = ^ . (3.4c) 

|Vc(r')|^x(r-r') 



V 



The next step in the coarse-graining procedure is the choice of the nonequihbrium en- 
semble. The natural choice is for it to be of the mixed ensemble due to the choice of the 
macroscopic variables 



41| . Since the hydrodynamic fields are simply mapped from the 
Cahn-Hilliard to the Doi-Ohta level, the appropriate probability distribution function is of 
the generalized microcanonical type. For the configurational variables, on the other hand, 
we choose a generalized canonical ensemble with the corresponding Lagrange multipliers Aq, 
Aq and A;. The total probability measure Py[x\ then takes the form 

Py[x\ = 5(np-p)5(ng-g)5(n, -e)p(Q,q)[c], (3.5a) 



P(Q,q) 



C 



N[y] 

where the normalization factor N[y] is 



exp 1^-^ rfVA(r) : n[c](r)j , (3.5b) 



N[y] 



j I?c^]l[x]exp j (fr\{r) : n[c](r)^ , (3.5c) 



and VLi[x\ = exp (^S^^\x)/kB). The projection operators Uq, Ilq, and 11; are, for simplicity, 
denoted by 11 in the above equations. The Lagrange multipliers Aq, Aq, and A;, denoted by 
A, are determined by the values of the slow variables y = (x) where the average is performed 
according to (13. ip . Their interpretation and the identification of the pro per values for the 
situation of interest is difficult and requires dynamic material information 48, 49|. However, 
for the results presented in this paper, the exact form of the nonequihbrium ensemble will 
not be needed. 
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B. Energy 



The energy of the coarse-grained Doi-Ohta level is obtained by averaging the energy of 
the more detailed Cahn-Hilliard level, 



j^dh jvxpy[x] l^dVx(r-r') (^^^ + 6 (r') + ^/.e |Vc(r')|' j . (3.6) 

Since we assume that the hydrodynamic fields p, g, and e are smooth on the more microscopic 
Cahn-Hilliard length scale, the energy expression takes the form 



Jv 

where in the last term on the right hand side of the equation, we have recognized the new 
variables Q(r) and /(r), determined by the mapping f l3.4al) and f l3.4cp . The last term on 
the right hand side of the above equation describes the energy density due to the presence 
of the interface which is proportional to the interfacial tension a. We therefore obtain the 
following well-known expression for the interfacial tension 

-(r) = (3.8) 
which here follows directly as the first result of the performed coarse-graining procedure. 



C. Entropy 

The coarse-grained entropy for the case of the generalized mixed ensemble Py[x\ is ob- 
tained from 

S^'\y) = J Pcp(Q,q)[c] [S^'\x) - A;b In p(Q,q)[c]] , (3.9) 

where S^^\x) is given by ( ]2.4bl) . The additional entropy in the coarse-grained expression, 
beside a simple ensemble average of the entropy from the more detailed level 1, is associated 
with the passage from the more microscopic configurational variable c, to the more macro- 
scopic variables Q, q, and /, while the hydrodynamic variables are taken to the coarser level 
without affecting the entropy. The additional entropy takes into account all the microstates 
with the composition c(r) consistent with the more coarse-grained state given with Q{r), 
q(r), and /(r). The functional derivative of the coarse-grained entropy is 

0,— -,AQ(r), Aq(r),Ai(r) . (3.10) 



5y \ dp ' ' T(r) 
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As discussed above, there exists a systematic procedure to obtain the Lagrange multiphers 



48| . However, their determination is a 



from the thermodynamically guided simulations 
difficult task and for the calculation of the time friction matrix elements in section UlI El we 
will assume A = 0. 



D. Poisson operator 



The general expression for the coarse-grained Poisson operator is given by 



L^'^ {y) 



'5 11 [x) 
6x 



6U {x 
6x 




(3.11) 



which for coarse graining from the Cahn-Hilliard to the Doi-Ohta level takes the form 



x/.V3/'.V.^54^4i'K^.)«^. (3.12) 

Jv Jv oxkirs) oxi{ri) 

where the mappings 11 are given by (13. 4p . and the elements of the Poisson operator L'-^^ 
by (12. 6p . To understand the above expression, note that the integrations over and r4 
come from the contraction of the operator L^^^ with 511/ 5x, which includes both matrix 
multiplication and integration over the position label. The obtained quantity must then be 
averaged in order to obtain the expression at the Doi-Ohta level. Therefore, the integrations 
over r[ and r'j come from the spatial smoothing of this quantity, since the ensemble averaging 
also implies smoothing in space, as noted in the previous subsection. 

From the elements of L^^^ in (12. 6 p and the mappings 11 in (13. 4p . we conclude that the 
Poisson operator L*^^) has the following form 

/ Li?} \ 



L(2)| 



ri,r2) 



.(2) .(2) .(2) .(2) .(2) .(2) 
l^trn J^gg J^e-f J^„D J^trn 





L 



-(2) 
(2) 



41 
i!? 

































(3.13) 
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In the analytical derivation of the coarse-grained Poisson operator, several approximations 
that are related to the difference in length scales (13. 2 p are used. The crossover in length scale 
from the more microscopic Cahn-Hilliard level to the more macroscopic Doi-Ohta level is 
performed through the smoothing function x(ri — r2). While the details of the derivation for 
some of the elements of L^^^ are given in the Appendix A, here we give their final expressions 



g(r2) 



dx dx 



dvi 



dx 



dv2 

dx 



'(2) 



d 



iri,r2j 



dr2a 
d 

dr2a 
d 

2/3 



2Pa/3(r2) - ^i^(r2)5a/3^ X 



dx 
dr 



2/3 



(3.14a) 
(3.14b) 
(3.14c) 
(3.14d) 



P,37(r2) + -P(r2)5/37 X 



dr 



Pa^ij2) + -P(r2)5a7 ) X 



2a ^ 



X 



d 



dr 



-Pa/3ir2j X 



L 



(2) 



.ri,r2j 
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^(ri) , ^ dx 



3 dr2p 



dl{Yi 



-X- 



(3.14e) 
(3.14f) 



P(ri)^"''''' (9r2a 3 dr2p drip 
where x = x{^i " ^2)- With this, we obtained all the elements of the coarse-grained Poisson 
operator since the elements L^gp-, Lgp, and -^g^'' arc obtained from the antisymmetry 
condition for L^'^\ We note the occurence of the smoothing function x(ri — r2) in the 
Poisson operator L^'^\ Indeed, when looking at the elements of L^^^ which involve only 
the hydrodynamic fields p, g, and e, we see that they differ from the appropriate elements 
of L^^^ in (12. 6p only in locality, i.e., instead of the Dirac delta function, we rather have a 
smoothing function x(ri — r2). While the original Cahn-Hilliard model is local in space, 
the coarse-grained model, instead, takes into account the whole volume element v{ri) of the 
smoothing function x(ri — r2). By assumption of the length scales comparison (13.21) . this 
smoothing function behaves simply as a delta function on the Doi-Ohta level due to the 
difference in length scales. Alternatively, the expressions for the elements of the Poisson 
matrix (I3.14dp - (13.14f]) could be also obtained based only on the mathematical character of 
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the variables P, p, and /. Similar to the original Doi-Ohta derivation, one would consider 
the transformation properties of the vector Vc and the mappings (13 ■4p . in order to infer the 
convective behavior. 

The Poisson operator L^^^ determines the convective behavior of the state variables, 
which can then be compared to the original Doi-Ohta model. The reversible time evolution 
equations for the variables {p, g, e, P, p, /} are obtained as 

%(ri) 

rev V 



dt 



ij 1^1, r2j 



%(r2) 



(3.15) 



under the previously discussed assumption that the smoothing function x(ri— r2) is acting as 

a Dirac delta function 5(ri— r2) at the Doi-Ohta level of description. Then, by transformation 

of the variables, the full set of the reversible time evolution equations for the set of state 

variables {p, g, e, Q, q, /} takes the form 

dp 
di 
dg 
dt 
de 
di 
dQ 
dt 
(9q 
'dt 



-V ■ (pv) , 






(3.16a) 


-V • (vg) - 


- Vp- V 


■r(^q-^gi^ -Qvr, 


(3.16b) 


-V ■ (ev) - 


P(V- v). 




(3.16c) 


-V ■ (gv) - 


- (Vv)^ : 


q+^Q(V-v), 


(3.16d) 


- (Vv) ■ q - 


-q.(Vv) 


^-^Q7-v- Vq+il(Vv) :q 





m 

dt 



+igi(v-v) + i (^q + igi ) ( q 

-V ■ (/v) + ^q : (Vv)^ + ^/ (V ■ v) . 



3^1 ) : (Vv) 



(3.16e) 
(3.16f) 



These equations correspond to the convective part of the time evolution equations of the 
Doi-Ohta model expressed in the GENERIC formalism Indeed, when the closure ap- 
proximation (11. 5p for the fourth moment of the interfacial normal vector nnnn is used, one 
obtains the exact equations as above (see Eqs. fll.4al) and (11. 6p ). Furthermore, since the Ja- 
cobi identity of the starting L^^^ operator is fulfilled and the closure approximation we used 
corresponds to the one made by Doi and Ohta 28|], and examined in 33|, we assume that 
the Jacobi identity of the derived Poisson operator L^"^^ is valid. However, this assumption 
might be questioned due to the approximations used in its derivation and which take into 
account different length scales (see Appendix A) . 
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E. Irreversible dynamics 



In this section we turn to the dynamic material properties which are contained in a 
friction matrix. There are two contributions to the coarse-grained friction operator, M^"^^ = 
M^^)' + M*^^)" . The first contribution, M^^^' , is obtained directly by averaging the elements 
of the friction matrix M^^\ in the same way as for the Poisson operator in (13.111) . This 
direct contribution describes dissipative effects that are carried on from a more microscopic 
level of description. A second contribution to the coarse-grained friction operator, M*^^)", 
results from the processes that are slower than the characteristic time scale of the Cahn- 
Hilliard level (given by the diffusion time ^d) but are fast compared to the processes at the 
Doi-Ohta level of description (with the time scale of the interface relaxation time Tdo,i,2)- 
This contribution presents an important feature of the coarse-graining procedure, since it 
captures the additional dissipation arising from the additional processes which can be treated 
as fluctuations on the time scale of the Doi-Ohta level. It can be evaluated from the Green- 
Kubo formula 



where Ts separates the times scales between the fast and slow variables, 11^ is the rapidly 
fluctuating part of the time derivative of the microscopic expressions for the slow variables 
y, and the average is over the atomistic trajectories consistent with the coarse-grained state 
y at t = and evolved according to the microscopic dynamics to the time t |41|. 

In order to identify the main dissipative contribution at the coarse-grained level and its 
origin, one has to understand all the processes occurring in a phase separation of a binary 
fluid under shear flow. In particular, through numerical simulations we tried to understand 
if there exist any processes in this system that can be considered as fast compared to the 
time scale of the Doi-Ohta level, and to estimate their contribution to dissipation at the 
Doi-Ohta level through the Green-Kubo part of the friction matrix M*^^-*", Eq. (I3.17p . We 
performed simulations of the Cahn-Hilliard equation fl2.9d|) for a binary mixture of volume 



Ts 



fraction = 0.3 which was subjected to a shear flow at time to = 200tD after the start of a 
phase separating process (when the well-defined droplets were formed). We have used the 
c^-structure for the entropy density s, as given in subsection III A[ We used a 1-dimensional 
equilibrium profile between the two coexisting bulk phases, with the composition field c{x) = 
sj a/6tanh(a;/2^), in order to define the length and time scale using ^ = \fK!^Ja and to = 
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100 200 300 400 % 2 4 6 

FIG. 1: (Color online) Time evolution of the coarse-grained Doi-Ohta variables for a binary mixture 
with (p = 0.3, subjected to the shear flow at time Iq = 200tD after the beginning of a phase 
separating process: interfacial width / (upper left), average amount of interface per unit volume Q 
(lower left), and shear viscosity, rj = —Qxy/j as a function of strain 7 = 7t (right). The results for 
the shear rates 7^0 = 0.05 and jto = 0.06 are presented. 

/Ma. Hydrodynamic interaction was not included and, therefore, coarsening of domains 
proceeded governed by the diffusion process, i.e., through the evaporation-condensation 
mechanism and the growth law L(t) ~ t^/^js], [sil, 32]. Once the well defined interfaces were 
formed, we looked for the fast processes permanently at work (like deformation of interface 
shapes through fast coagulation and break-up, etc.), which would give rise to new irreversible 
dynamics, going beyond the diffusion and hydrodynamic ones. 

The time evolution of the coarse-grained Doi-Ohta variables in FigH] is computed by a 
simple use of mappings (13.41) . The shear viscosity, which arises from excess shear stress due 
to the domain interfaces, is related to the anisotropy element q^y through rj = —qxy/j- From 
the time evolution of the average interfacial width /, one can see that this new structural 
variable relaxes very quickly and is afterwards only affected by the flow through the convec- 
tive behavior, Eq.( l3.16fj) . This numerical analysis shows that the coarse-grained Doi-Ohta 
variables fluctuate in time (with the time scale of ~ I/7) around their average values due 
to the competition between the flow field and the ordering mechanisms, see Fig{Tl This can 
be also seen in the time evolution of the morphology. The snapshots of the morphology, 
presented in Fig|2l show that at first domains are elongated due to the shear flow, which 
tends to orient them towards the shear direction. The elongation is then followed by coagu- 
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7 = 7 = 3 7 = 6 7 = 9 7 = 12 



FIG. 2: Configurations of a binary fiuid with volume fraction (p = 0.3, simulated on a 128 x 128 
lattice during phase separation under shear flow at different values of strain 7 = 'yt. Shear flow 
was applied after initial time of to = 200tDj with the shear rate = 0.05. The snapshots were 
taken until the maximum strain 7 = 12. The gray scale is used in order to capture the processes 
of interface deformation. 



ets, 



lation, break-up, shape relaxation, as well as the evaporation-condensation of the drop 
which are exactly the mechanisms of interfacial relaxation proposed by Doi and Ohta 
and Lee and Park 3J| . Similar interfacial relaxation processes have been identified in other 



28|, 



works which report simulations of the Cahn-Hilliard model with and without hydrodynamic 



mteraction 
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2J]. The morphology time evolution shows no new processes that 
are fast compared to the Doi-Ohta level of description, and which would emerge as dissi- 
pation on the coarse-grained level through the Green-Kubo formula fl3.17p . Moreover, the 
time evolution curves of the Doi-Ohta variables in Fig. [1] are smooth without the presence 
of thermal fluctuations (which were not incorporated into our analysis). Since the above 
described mechanisms of interfaces relaxation are also the ones considered for the analysis 
of different growth regimes in the late stage kinetics of phase separation (see Section [T]), we 
conclude that the main dissipative contribution at the coarse-grained level arises through 
the averaging of the already present dissipative terms at the Cahn-Hilliard level, M*^^^'. 

The relaxation of the Doi-Ohta variable y = {p, g, e, P, p, 1} could be expressed from the 
irreversible part of the GENERIC time evolution equation (12. 2p as 



%(r 



dt 



c^V.<(r„r.)fi^. (3.18) 
^ %(r2) 



In order to derive the coarse-grained friction matrix M^"^^' , we note that the Cahn-Hilliard 
friction matrix M^^^ consists of a hydrodynamic and a diffusive part, M^^^ = M^^^'^^'^ + 
given in subsection III A[ Therefore M*^^-*' will also comprise the hydrodynamic 
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part and a part arising from the diffusion. Since the hydrodynamic variables are simply 
taken over from the Cahn-Hilliard level, the appropriate hydrodynamic elements of the new 
friction matrix are the same as in fl2.7p . The diffusive part, on the other hand, is related to 
the configurational variable c, and the ensemble averaging must be performed in a similar 
way as for the Poisson operator in Section IIIIDI For the analytical derivation of the coarse- 
grained friction matrix M'^^^'''^'^, in analogy to Eq. (13. lip , one starts from 

M(^)''- (,) = / ' . M(^)-- (X) . \ , (3.19) 

\ I y 

and then employs similar procedure and approximations as in the derivation of the elements 
of the Poisson matrix L*^^) (see Appendix A). For example, the relaxation of the Doi-Ohta 
variable Q could be expressed as an irreversible part of the GENERIC time evolution equa- 
tion fIXISl) as 

^(ri) / dh^ Mj,f ■'^^Vi, 1-2) (3.20) 



ag(ri) 



where we have used the approximation {dQijc) / dt)\^^^^^ ~ /(r) {dPij) / dt)\^^^^^, since the 
interfacial width / does not change much, compared to the other two configurational vari- 
ables. Using the entropy functional derivative f l3.10p with A = 0, results in the matrix 

Ciy dif 

element Mp^ ' being the only relevant term for the variable P which is left after the av- 
eraging f l3.19p . Then, putting the equation f l3.20p in the Doi-Ohta equation f ll.7ap for the 
relaxation of Q in a homogeneous case gives the following expression for the relaxation time 



1 _ 2keMI 

Tdo,l VQ 



^cp(Q,q)[c] / d^r 
Jv 



dr \ dr"^ 



where V is the volume of the system. Details of the calculation are given in Appendix B. 
Similar expression could be obtained for the relaxation time rdo,2- The above formula, gives 
the expression for the Doi-Ohta relaxation time t^o i calculated from the transformation 
of the element of the friction matrix M^^^'^^^ containing the relaxation parameters of the 
Cahn-Hilliard level at any time t. This is the result of the presented systematic coarse- 
graining procedure developed within the GENERIC formalism which shows in which way 
the Doi-Ohta relaxation times can be obtained from the more detailed level of description. 
However, the above relation is difficult to express analytically in form of the variables of 
the Doi-Ohta level, due to the diffusive terms coming from the Cahn-Hilliard level. On the 
other hand, the above expression could, in principle, be calculated numerically using short 
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time simulations at the Cahn-Hilliard level. However, the problems concerning the third 
order derivatives of the composition field c make this task rather complicated, as well as the 
need for determination of the Lagrange multipliers in order to employ the full GENERIC 
procedure developed in this paper. 

IV. CONCLUSIONS AND PERSPECTIVES 

In this paper, we have developed a systematic, thermodynamically guided method which 
establishes the coarse-grained Doi-Ohta model, used for rheological behavior, from the more 
detailed Cahn-Hilliard model of a phase separating binary fiuid. The main contributions 
of the present work can be summarized as: (1) the introduction of the average interfacial 
width I as an additional structural variable; (2) the derivation of the coarse-grained Poisson 
operator L^'^\ Eg. (13. lip , from its finer level analogue; (3) no new dissipative processes arise 
during level jumping apart from the known hydrodynamic and diffusive dissipation; (4) the 
derivation of an expression fl3.2ip for the Doi-Ohta relaxation times as a function of the 
elements of the friction matrix describing relaxation at the finer Cahn-Hilliard level. 

The crucial step in the coarse-graining procedure was to identify the relevant coarse- 
grained variables and find the appropriate mapping which expresses them in terms of the 
more microscopic variables. In order to capture the physics of the Doi-Ohta level, we in- 
troduced the interfacial width l{r,t) as an additional variable in the new model. In that 
way we could account for the stretching of the interface under fiow and derive analytically 
the reversible (convective) behavior of the variables Q{r,t) and q(r,t), which recovers the 
already established phenomenological Doi-Ohta model. Introduction of the interfacial width 
I and its time evolution equation, as an addition to Q and q of the original Doi-Ohta model 
with zero-width interface, offers new possibilities for modeling of the rheological behavior 
of multiphase flows. In addition, the expression for the interfacial tension f l3.8p in terms 
of the Cahn-Hilliard parameters follows as the direct result of the developed systematic 
coarse-graining procedure. 

Considering the irreversible dynamics on the Doi-Ohta level, it has been shown that the 
dissipative processes at this coarse-grained level are carried on from the more microscopic 
Cahn-Hilliard level. Although their analytical derivation is too complex and rich in structure, 
the way to connect the interface relaxation times of the Doi-Ohta model and the underlying 
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morphology dynamics simulated at the Cahn-Hilliard level is established. Analysis of the 
numerical investigation of phase separation under shear in the diffusive regime revealed no 
new physical process that occur at the time scale which is slow compared to the Cahn- 
Hilliard level and fast from the perspective of the Doi-Ohta level. Therefore, there are no 
new dissipative effects arising in the performed coarse-graining step through the Green- 
Kubo type formula, i.e., no new physics appear. That leads us to the conclusion that all 
the dissipation has been introduced into the model by coarse graining from the reversible 
atomistic level to the Cahn-Hilliard level. This coarse-graining step has been done in the 
derivation of the hydrodynamic equations of a phase separating fluid mixture from the 
underlying microscopic dynamics 46|. Furthermore, Espahol and Vazquez showed in 50 1 
that all dissipation is introduced already at the very first coarse-graining step, going from 
the atomistic to the Fokker-Planck level. From there then follows the conclusion that the 
Cahn-Hilliard level is not so fundamental considering that the same dissipation mechanisms 
occur even at the finer levels, and are afterwards transmitted to an even more coarse-grained 
Doi-Ohta level, with no new dissipative processes arising. 

The presented results could be used to deal with several interesting open problems. First, 
possible extension of the procedure to more complicated phase separating systems with the 
addition of surfactants or block copolymers would be very interesting for wider use. Second, 
the nature of the dissipative processes in a phase separating binary mixture is still far 
from clear. While during the diffusive regime, in the absence of thermal noise, we do not 
recognize any further fast processes that could be considered as fluctuations at the Doi- 
Ohta level, this becomes far from evident during the late stages of coarsening when inertia 
becomes important. The nature of the domain coarsening in the inertial hydrodynamics 
regime is still unknown, as well as the asymptotic behavior of a phase separating system 
under shear. Furthermore, recent studies showed that a nonequilibrium steady state can 
be reached at high Reynolds numbers (low shear rates) only in the mixed viscous-inertial 
regime, i.e., with any non-zero amount of inertia 2ll-l24|. The fact that no steady state could 
be formed without presence of inertia, no matter how small, suggests that inertia plays the 
role of a singular perturbation in this problem The formed steady state, arising from a 
cyclic occurrence of elongation of domains followed by their break up, coagulation and shape 
relaxation, is characterized by irregular fluctuations around the attained steady state values. 
It would be interesting to understand the nonequilibrium steady state and the origin of these 
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fluctuations, using the presented coarse-graining procedure. However, the phenomenological 
Doi-Ohta theory assumes low Reynolds number, i.e., neglects inertia. Deeper insight into the 
mixed viscous- inert ial regime is therefore crucial in order to develop the connection between 
the two levels and to understand the role of inertia in this problem. 



Appendix A 

(2) (2) 

We present a derivation of the elements Lpg(ri, and Ly (ri, of the coarse-grained 
Poisson operator L^^-* in detail. For the element Lpg(ri,r2) the general expression fl3.12p 
gives 

4'g(ri,r2) = jvxpy[x\ j dh[ j X(ri - r'l) x(r2 - r'2) 

Jv Jv ^c(r3) ^ 5g(r4) 

which after substitution of the appropriate functional derivatives. Leg element, and the 
integration over and becomes 

Lg^(ri,r2) = !vcp(Q,^)[c] j dh[ [ d'r', xi^, - r[) xir 2 - r',) 
J Jv Jv 



dc d 



For further calculations, we use the difference between the Cahn-Hilliard and Doi-Ohta 



length scales (13. 2p . so that fy d^r[ xi'^i " — ~ x(ri — r'g) (see (45|). Then the last 

equation, after integration over t[, becomes 

4'l(ri, r2) = / Pcp(Q,q)[c] 1^ rfM X{r. " ^2)^ " (2 ^ x{r^ - r'2)) . (A3) 

After integration by parts, in which surface terms vanish, using approximation x(ri — 
^2)x{^2 — ^2) ~ xi^i ~ r2)x(r2 — i'2) valid due to the fact that the smoothing length is much 
smaller than the Doi-Ohta length scale, equation (13. 2p . and the identity (9x(r — r')/dr' = 
— 5x(r — r')/dr, we come to the expression 



/ d^r'2 |Vc(r^)|'x(ri-r2)x(r2-r'2) 
Jv 



d 



dr 



2/3 



j^dh'2 |Vc(r^)|'x(ri-r2)x(r2-r;,)| 



(A4) 
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Finally, with the help of the mappings (13. 4p . we can recognize the appropriate terms in the 
last equation, so that their ensemble average gives 



d 



dr- 



2a 



PM ^^t'"''' . (A5) 



dr- 



2/3 



The Poisson operator element L\'^{ri,r2) is obtained under similar approximations as 
above, but also the following approximations 

[ d'r[ ^c]{r[) x(ri - v[) xir[ - v',) ^ ^c]{r^) xir^ - r^), (A6) 
Jv 

which holds under the assumption of the difference in length scales ( 13. 2p . and the assumption 

((np[c](r))"n,[c](r)) ^ ((np[c](r)))" (n,[c](r)) , 
for n = ±1. For the element Ljg (ri, r2), we then obtain 



(A7) 



4?«(ri,r2) 



d^r' J^^|Vc(r')r'x(ri-r') 



dx{ri - ra) 



P{ri) \JV 

dKri) 2 9x(ri-r2) 
— ^ X ri - r2 - -/ ri 

o ^ri) / 9x(ri-r2) 



9r 



2a 



(A8) 



In order to express the previous equation in terms of the coarse-grained Doi-Ohta variables 
{P, p, /}, we must use the following closure approximation for the first term on the right 
hand side 



£r' ^^|Vc(r')r'x(ri-r') 



V 



C 



dc dc 



V 



-X[ri - r 



/ rfV |Vc(r')|x(ri -r'^ 
Jv 



(A9) 



where C is chosen in such a way that the previous expression is valid when its trace is taken. 
By putting a = P into the above equation, we find 



C 



d-^r' |Vc(r')rx(ri -r'^ 



(AlO) 



Using the above expressions, element Lj^ (ri,r2) takes the final form 



^Syri,r2) 



/(ri) 9x(ri-r2) 1 dxir^ - ^2) 

dl{r 
dr 



— X(ri - r2j 



1/3 



(All) 
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Appendix B 



We present the derivation of the expression fl3.2ip for the time relaxation rdo,i- Substitu- 
tion of the equation fl3.20p for the irreversible time evolution of the configurational variable 
Q into the Doi-Ohta relaxation equation (II .yap gives 

Tdo,! <5(ri) Jv T{r2) 

Under the assumption of a homogeneous system, integration over the whole system size of 
the above expression gives 



1 / 



rfVi / dh^MP^'-^^vuv^). (B2) 
V Jv 



Tdci QVT 

The general expression (13.190 for the coarse-grained friction matrix M'^^^'''^'^ gives 
''^'^(ri,r2) = Jvxpy[x] J d^r[ j d^r'^xiji - ^i) xij2 - r'^) 

Jv Jv c'c(r3) ()e(r4) 

which after substitution of the appropriate functional derivatives, M^^^''^^^ element, and the 
integration over r^, r[, and r'g becomes 



where we have used the assumption of a homogeneous (M = const.) and isothermal system, 
as well as the approximations based on a difference between the Cahn-Hilliard and Doi- 



Ohta length scales, as in Appendix A. Substitution of the friction matrix element (lB4p into 
equation flB2p . gives after using the normalization condition jyd^r' x(r — r') = 1, the final 
formula for the Doi-Ohta time relaxation 



2keMI 



Tdo,l VQ 



^cp(Q,q)[c] / d^r 
Jv 



d ( d'^c(r) 



dr V dr"^ 



(B5) 
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